AFRL-OSR-VA-TR-201 4-0232 


Design  Optimizations  Simulation  of  Wave  Propagation  in  Metamaterials 


Robert  Freund 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 


09/24/2014 
Final  Report 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Air  Force  Research  Laboratory 
AF  Office  Of  Scientific  Research  (AFOSR)/  RTA 
Arlington,  Virginia  22203 
Air  Force  Materiel  Command 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

14-09-2014  Final 

3.  DATES  COVERED  (From  -  To) 
15/06/2011-14/06/2014 

4.  TITLE  AND  SUBTITLE 

Design  Optimization  and  Simulation  of  Wave  Propagation  In 

Metamaterials 

5a.  CONTRACT  NUMBER 

FA9550-11-1-0141 

5b.  GRANT  NUMBER 

FA9550-1 1  -1-0141 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Freund,  Robert 

Peraire,  Jaime 

Nguyen,  Cuong 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Massachusetts  Institute  of  Technology 

77  Massachusetts  Avenue 

Cambridge,  MA  02139 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Science  and  Research 

875  Randolph  Street 

Suite  325  Room  3112 

Arlington,  VA  22203 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFOSR 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Metametarials  are  scientifically  engineered  materials  that  are  designed  to  interact  with  and  control  mechanical  and 
electromagnetic  waves  in  ways  that  cannot  be  achieved  with  conventional  materials.  For  instance,  metamaterials  can  be 
designed  to  bend  electromagnetic  waves  around  an  object  so  that  the  object  appears  invisible  to  surrounding  observers,  focus 
light  to  create  a  subwavelength  image  of  a  source,  and  mitigate  blast  waves  to  protect  structures  and  humans  from  explosion. 
Among  the  many  research  fronts  we  have  been  working  on  under  this  grant,  one  highlight  is  the  introduction  of  a  fundamentally 
new  modeling  paradigm  that  incorporates  fabrication  issues  in  the  optimal  design  of  metamaterials.  Another  highlight  of  our 
work  under  this  grant  is  our  successful  development  of  full  3D  topology  optimization  (in  which  every  voxel  of  the  unit  cell  is  a 
degree  of  freedom)  of  photonic-crystal  structures  in  order  to  find  optimal  omnidirectional  band  gaps  for  various  symmetry 
groups,  including  fee  (including  diamond),  bcc,  and  simple-cubic  lattices. 

15.  SUBJECT  TERMS 

wave  propagation,  metamaterials,  periodic  media,  band-gap  optimization,  semidefinite  programming,  photonic  crystal  design,  phi 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 


Prescribed  by  ANSI  Std.  Z39.18 
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Design  Optimization  and  Simulation  of  Wave  Propagation  in  Metamaterials 


Metametarials  are  scientifically  engineered  materials  that  are  designed  to  interact  with  and 
control  mechanical  and  electromagnetic  waves  in  ways  that  cannot  be  achieved  with 
conventional  materials.  For  instance,  metamaterials  can  be  designed  to  bend  electromagnetic 
waves  around  an  object  so  that  the  object  appears  invisible  to  surrounding  observers,  focus  light 
to  create  a  subwavelength  image  of  a  source,  and  mitigate  blast  waves  to  protect  structures  and 
humans  from  explosion. 

Among  the  many  research  fronts  we  have  been  working  on  under  this  grant,  one  highlight  is  the 
introduction  of  a  fundamentally  new  modeling  paradigm  that  incorporates  fabrication  issues  in 
the  optimal  design  of  metamaterials.  This  new  modeling  method  is  motivated  by  the  observation 
that  the  computed  optimal  solution  of  an  optimization  problem  often  cannot  be  implemented 
directly,  irrespective  of  data  accuracy,  due  to  either  (i)  technological  limitations  (such  as  physical 
tolerances  of  machines  or  processes),  (ii)  the  deliberate  simplification  of  a  model  to  keep  it 
tractable  (by  ignoring  certain  types  of  constraints  that  pose  computational  difficulties),  and/or 
(iii)  human  factors  (getting  people  to  “do”  the  optimal  solution).  We  have  developed  a  modeling 
paradigm  called  “fabrication-adaptive  optimization”  for  treating  issues  of 
implementation/fabrication.  We  develop  computationally-focused  theory  and  algorithms,  and 
we  present  computational  results  for  incorporating  considerations  of  implementation/fabrication 
into  constrained  optimization  problems  that  arise  in  photonic  crystal  design.  The  fabrication- 
adaptive  optimization  framework  stems  from  the  robust  regularization  of  a  function  (see  A. 

Lewis  (2002)  “Robust  regularization,”  Technical  Report,  Simon  Fraser  University,  Burnaby, 
British  Columbia,  Canada).  When  the  feasible  region  is  not  a  normed  space  (as  typically 
encountered  in  application  settings),  the  fabrication-adaptive  optimization  framework  typically 
yields  a  non-convex  optimization  problem.  (In  the  special  case  where  the  feasible  region  is  a 
finite-dimensional  normed  space,  fabrication-adaptive  optimization  can  be  re-cast  as  an  instance 
of  modem  robust  optimization.)  We  study  a  variety  of  problems  with  special  structures  on 
functions,  feasible  regions,  and  norms,  for  which  computation  is  tractable,  and  develop  an 
algorithmic  scheme  for  solving  these  problems  in  spite  of  the  challenges  of  non-convexity.  We 
apply  our  methodology  to  compute  fabrication-adaptive  designs  of  two-dimensional  photonic 
crystals  with  a  variety  of  prescribed  features.  The  fabrication-adaptive  methodology  together 
with  the  convex  optimization  fonnulations  we  proposed  before  integrate  into  a  systematic  design 
framework  for  the  metamaterials  applications  that  are  of  direct  and  practical  interest  to  the  Air 
Force. 

Another  highlight  of  our  work  under  this  grant  is  our  successful  development  of  full  3D  topology 
optimization  (in  which  every  voxel  of  the  unit  cell  is  a  degree  of  freedom)  of  photonic-crystal 
structures  in  order  to  find  optimal  omnidirectional  band  gaps  for  various  symmetry  groups, 
including  fee  (including  diamond),  bcc,  and  simple-cubic  lattices.  Even  without  imposing  the 
constraints  of  any  fabrication  process,  the  resulting  optimal  gaps  are  only  slightly  larger  than 
previous  hand  designs,  suggesting  that  current  photonic  crystals  are  nearly  optimal  in  this 
respect.  However,  optimization  can  discover  new  structures,  e.g.  a  new  fee  structure  with  the 
same  symmetry  but  slightly  larger  gap  than  the  well-known  inverse  opal,  which  may  offer  new 
degrees  of  freedom  to  future  fabrication  technologies.  Furthermore,  our  band-gap  optimization  is 


an  illustration  of  a  computational  approach  to  3D  dispersion  engineering  which  is  applicable  to 
many  other  problems  in  optics,  based  on  a  novel  semi  definite  program  formulation  for  non- 
convex  eigenvalue  optimization  combined  with  other  techniques  such  as  a  simple  approach  to 
impose  symmetry  constraints.  We  also  demonstrate  a  technique  for  robust  topology  optimization, 
in  which  some  uncertainty  is  included  in  each  voxel  and  we  optimize  the  worst-case  gap,  and  we 
show  that  the  resulting  band  gaps  have  increased  robustness  to  systematic  fabrication  errors. 
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